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Abstract. Group testing algorithms are very useful tools for DNA 
library screening. Building on recent work by Levenshtein (2003) and 
Tonchev (2008), we construct in this paper new infinite classes of combi- 
natorial structures, the existence of which are essential for attaining the 
minimum number of individual tests at the second stage of a two-stage 
disjunctive testing algorithm. 

1 Introduction 

With the completion of genome sequencing projects such as the Human Genome 
Project, efficient screening of DNA clones in very large genome sequence data- 
bases has become an important issue pertaining to the study of gene functions. 
Very useful tools for DNA library screening are group testing algorithms. The 
general group testing problem (cf. [9]) can be basically stated as follows: a large 
population X of v items that contains a small set of defective (or positive) items 
shall be tested in order to identify the defective items efficiently. For this, the 
items are pooled together for testing. The group test reports "yes" if for a subset 
SCI one or more defective items have been found, and "no" otherwise. Using 
a number of group tests, the task of determining which items are defective shall 
be accomplished. Various objectives could be considered for group testing, e.g., 
minimizing the number of group tests, limiting the number of pools or pool sizes, 
or tolerating a few errors. In what follows, we will focus on the first issue. 

Of particular practical importance in DNA library screening are one- or two- 
stage group testing procedures (cf. [18, p. 371]): 

"[...] The technicians who implement the pooling strategies generally dislike 
even the 3-stage strategies that are often used. Thus the most commonly used 
strategies for pooling libraries of clones rely on a fixed but reasonably small 
set on non-singleton pools. The pools are either tested all at once or in a small 
number of stages (usually at most 2) where the previous stage determines 
which pools to test in the next stage. The potential positives are then inferred 
and confirmed by testing of individual clones [...]." 

* The author gratefully acknowledges support of his work by the Deutsche Forschungs- 
gemeinschaft (DFG) via a Heisenberg grant (Hu954/4) and a Heinz Maier-Leibnitz 
Prize grant (Hu954/5). 
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Disjunctive testing relies on Boolean operations. It aims to find the set of de- 
fective items by reconstructing its binary (0, l)-incidence vector x = (sci, . . . , x v ), 
where Xi — 1 if the ith item is defective, and x% — otherwise. Levenshtein [20] 
(cf. also [25]) has employed a two-stage disjunctive testing algorithm in order 
to reconstruct the vector x: At Stage 1, disjunctive tests are conducted which 
are determined by the rows of a binary matrix that is comparable to a parity- 
check matrix of a binary linear code. After determining what items are positive, 
negative or unresolved, individual tests are performed at Stage 2 in order to 
determine which of the remaining unresolved items are positive or negative. 

Particularly important with respect to the research objectives in this paper, 
Levenshtein derived a combinatorial lower bound on the minimum number of 
individual tests at Stage 2. He showed that this bound is met with equality 
if and only if a Steiner i-design exists which has the additional property that 
the blocks have two sizes differing by one (i.e., k and k + 1). Relying on this 
result, Tonchev [25] gave a straightforward construction method for such designs: 
Suppose that T> = (X,B) is a Steiner t-(v,k,l) design that contains a Steiner 
(t - l)-(v, k, 1) subdesign V = (X,B'), where B' C B. Then, the blocks of V, 
each extended with one new point x $ X, together with the blocks of V that do 
not belong to V , form a Steiner t-(v + 1, {k, k + 1}, 1) design. Relying on specific 
balanced incomplete block designs (BIBDs), he constructed two infinite classes 
of such designs: A Steiner 2-(q e + 1, {q,q + 1}, 1) design exists for every prime 
power q and every positive integer e > 2, and a Steiner 2-(6et + 4, {3, 4}, 1) design 
for every positive integer a, based on resolvable BIBDs from affine geometries 
and Kirkman triple systems, respectively. Moreover, he constructed two infinite 
classes derived from certain Steiner quadruple systems. 

In this paper, we build on the work by Levenshtein and Tonchev and con- 
struct several further infinite classes of Steiner designs with the desired additional 
property. Our constructions involve, inter alia, resolvable BIBDs and cyclically 
resolvable BIBDs. As a result, we obtain efficient two-stage disjunctive group 
testing algorithms suited for DNA library screening. 

The paper is organized as follows: Combinatorial tools and structures which 
are important for our further purposes are provided in Section 2. Section 3 
presents a short overview of the previous combinatorial approaches and construc- 
tions by Levenshtein and Tonchev. Section 4 is devoted to our new combinatorial 
constructions. The paper concludes in Section 5. 

2 Combinatorial Structures and Tools 

Let A be a set of v elements and B a collection of fc-subsets of X . The elements 
of X and B are called points and blocks, respectively. An ordered pair V — (A, B) 
is defined to be a t-(v, k, A) design if each t-subset of X is contained in exactly 
A blocks. For historical reasons, a t-(v, k, A) design with A = 1 is called a Steiner 
t-design or a Steiner system. Well-known examples are Steiner triple systems (t = 
2, k = 3) and Steiner quadruple systems (t = 3, k — 4). A 2-design is commonly 
called a balanced incomplete block design, and denoted by BIBD(i>, k, A). It can 
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be easily seen that in a t-(v, k, A) design each point is contained in the same 
number r of blocks, and for the total number b of blocks, the parameters of a 
t-(v, k, A) design satisfy the relations 



bk = vr and r{k - I) = \j—^-{v - I) for t > 2. 



Example 1. Take as point-set 

X = {1,2, 3, 4, 5, 6, 7, 8, 9} 

and as block-set 

B = {{1,2, 3}, {4, 5,6}, {7, 8,9}, {1,4, 7}, {2, 5, 8}, {3, 6, 9}, 

{1,5, 9}, {2,6, 7}, {3, 4, 8}, {1,6,8}, {2, 4, 9}, {3, 5, 7}}. 

This gives a BIBD(9,3,1), i.e. the unique affine plane of order 3. It can be 
constructed as illustrated in Figure 1. 




Fig. 1. A BIBD(9, 3, 1). 



In this paper, we primarily focus on BIBDs. Let (X, B) be a BIBD(t>, k, A), 
and let a be a permutation on X. For a block B = {6i, . . . , bk} £ B, define B a := 
{b°, . . . , 6£}. If B a := {B a : B E B} = B, then a is called an automorphism of 
(X,B). If there exists an automorphism a of order v, then the BIBD is called 
cyclic. In this case, the point-set X can be identified with Z„, the set of integers 
modulo v, and a can be represented by a : i -> i + 1 (mod u). 

For a block i? = {&i, . . . , &&} in a cyclic BIBD(u, fc, A), the set B + i := {&i +i 
(mod v), . . . ,bh + i (mod u)} for i g is called a translate of £?, and the set 
of all distinct translates of B is called the orbit containing B. If the length of 
an orbit is v, then the orbit is said to be full, otherwise short. A block chosen 
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arbitrarily from an orbit is called a base block (or starter block). If k divides v, 
then the orbit containing the block 



is called a regular short orbit. For a cyclic BIBD(w, k, 1) to exist, a necessary 
condition is v = 1 or k (mod k(k — I)). When v = 1 (mod fc(fc — I)) all orbits 
are full, whereas if v = k (mod k{k — 1)) one orbit is the regular short orbit and 
the remaining orbits are full. 

A BIBD is said to be resolvable, and denoted by RBIBD(u, k, A), if the block- 
set B can be partitioned into classes IZi , . . . , lZ r such that every point of X is 
contained in exactly one block of each class. The classes IZi are called resolution 
(or parallel) classes. A simple example is as follows. 

Example 2. An RBIBD(9,3, 1). Each row is a resolution class. 



Generally, an RBIBD(/c 2 , k, 1) is equivalent to an afhne plane of order k. An 
RBIBD(w,3,l) is called a Kirkman triple system. Necessary conditions for the 
existence of an RBIBD(w, k, A) are X(v — 1) = (mod (k — 1)) and v = (mod 



If lit is a resolution class, define K? := {B a : B e Hi}. An RBIBD is 
called cyclically resolvable, and denoted by CRBIBD(?;, k, A), if it has a non- 
trivial automorphism a of order v that preserves its resolution {IZi, . . . TZ r }, i.e., 
{Til, ...K?} = {Tlx, . . . TZ r } holds. An example is as follows (cf. [12]). 

Example 3. A CRBIBD(21, 3, 1) is given in Table 1. The base blocks are {1, 4, 16}, 
{19, 20, 3}, {1, 11, 19}, and {0, 7, 14}. There are three full orbits and one regu- 
lar short orbit. Each row is a resolution class. One orbit of resolution classes is 
{Tvlo, ■ • ■ ,Hq}, and another orbit is {1Zq,1Z[ 7 1Z' 2 }. 

Mishima and Jimbo [21] classified CRBIBD(w, k, l)s into three types, ac- 
cording to their relation with cyclic quasiframes, cyclic semiframes, or cyclically 
resolvable group divisible designs. They can only exist when v = k (mod k(k— 1)). 

In a cyclic BIBD(i>, k, 1), we can define a multiset AB := {hi — bj : i,j = 
1, . . . , k; i j} for a base block B = {b\, . . . , &*,}. Let {Bi} i& j, for some index 
set /, be all the base blocks of full orbits. If v = 1 (mod k(k — 1)), then clearly 




Hi 



{1,2,3} {4, 5, 6} {7, 8, 9} 
{1,4,7} {2,5,8} {3,6,9} 
{1,5,9} {2,6,7} {3,4,8} 
{1,6,8} {2,4,9} {3,5,7} 



|J ABi = Z v \ {0}. 



The family of base blocks {Bi}i<=i is then called a (cyclic) difference family in 
Z„, denoted by CDF(w, k, 1). 



Table 1. Example of a CRBIBD(21, 3, 1). 
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TZ 
Hi 

n 2 

lis 
TZ 5 



{1,4, 16} {8,11,2} {15,18,9} 
{2,5,17} {9, 12,3} {16, 19,10} 
{3,6, 18} {10, 13,4} {17,20,11} 
{4,7, 19} {11,14,5} {18,0,12} 
{5,8,20} {12,15,6} {19,1,13} 
{6,9,0} {13,16,7} {20,2,14} 
{7, 10,1} {14,17,8} {0,3,15} 



{19, 20, 3} {5, 6, 10} {12, 13, 17} {0, 7, 14} 

{20,0,4} {6,7,11} {13,14,18} {1,8,15} 

{0,1,5} {7,8,12} {14,15,19} {2,9,16} 

{1, 2, 6} {8, 9, 13} {15, 16, 20} {3, 10, 17} 

{2,3,7} {9,10,14} {16,17,0} {4,11,18} 

{3, 4, 8} {10, 11, 15} {17, 18, 1} {5, 12, 19} 

{4, 5, 9} {11, 12, 16} {18, 19, 2} {6, 13, 20} 



n' 2 



{1,11,9} {4, 14,12} {7,17,15} 
{2,12,10} {5,15,13} {8,18, 16} 
{3,13,11} {6,16,14} {9,19,17} 



{10, 20, 18} {13, 2, 0} {16, 5, 3} {19, 8, 6} 
{11,0,19} {14,3,1} {17,6,4} {20,9,7} 
{12,1,20} {15,4,2} {18,7,5} {0,10,8} 



Let k be an odd positive integer, and p = 1 (mod k(k — 1)) a prime. A 
CDP(p, k, 1) is said to be radical, and denoted by RDF(p, k, 1), if each base 
block is a coset of the fc-th roots of unity in Z p (cf. [3]). A link to CRBIBDs has 
been established by Genma, Mishima and Jimbo [12] as follows. 

Theorem 1. If there is an RDF(p, k, 1) with p a prime and k odd, then there 
exists a CRBIBD{pk,k,l). 

The notion of resolvability holds in the same way for t-(v, k, A) designs with 
t > 2. Moreover, a Steiner quadruple system 3-(u, 4, 1) is called 2-resolvable if 
its block-set can be partitioned into disjoint Steiner 2-(v, 4, 1) designs. 

For encyclopedic references on combinatorial designs, we refer the reader 
to [2,8]. A comprehensive book on RBIBDs and related designs is [11]. Highly 
regular designs are treated in the monograph [15]. A recent survey on various 
connections between error-correcting codes and algebraic combinatorics is given 
in [16]. For an overview of numerous applications of combinatorial designs in 
computer and communication sciences, see, e. g., [6,7,17]. 

3 Combinatorial Approaches and Constructions by 
Levenshtein and Tonchev 

This section presents a short overview of recent combinatorial approaches and 
constructions by Levenshtein [20] and Tonchev [25]. 

3.1 Two-Stage Disjunctive Group Testing Algorithm 

Disjunctive group testing relies on Boolean operations in order to solve the 
problem of reconstructing an unknown binary vector x of length v using the 
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pool testing procedure [9]. Particularly, Levenshtein [20] (see also [25]) has 
employed a two-stage disjunctive testing algorithm to reconstruct the vector 
x = (x\, ...,x v ): At Stage 1, disjunctive tests are conducted which are deter- 
mined by the rows of a binary u x v matrix H = (hi t j) that is comparable to 
a parity-check matrix of a binary linear code. A syndrome s = (si, . . . , s u ) is 
calculated, where Sj is defined by 

V 

j'=i 

where V and & denote the logical operations of disjunction and conjunction. 
The system of u logical equations with v Boolean variables for reconstructing 
the vector x = (x%, . . . , x v ) does not have a unique solution in general. After 
determining what items are positive, negative or unresolved, individual tests are 
performed at Stage 2 in order to determine which of the remaining unresolved 
items are positive or negative. 

3.2 Minimum Number of Individual Tests 

Let X(v) be the set of all 2" subsets of the set X = {1, 2, . . . , v} and X t (v) — 
{x € X(v) : \x\ = t}. For a fixed t (1 < t < v) consider a covering operator 
F : X t (v) -t X(v) such that x C F(x) for any x € X t (v). Define 

V = {F(x) : x g X t (v)}. 

For any T, 1 < T < (Tj, consider the decreasing continuous function gt(T) = 
k + ^7-^(1 — a) where k and a are uniquely determined by the conditions 
T( t ) = a(^), fee {t,...,v}, and 1 — < a < 1. Using averaging and linear 
programing, Levenshtein [20] proved the following inequality: 

Theorem 2 (Levenshtein, 2003). 

£ |F(x)|>. 9t (P|), 

and t/ie bound is met with equality if and only ifT> is a Steiner t-(v, {k, k + 1}, 1) 
design. 

One of the main motivations for the above result is to minimize the number 
of individual tests at the second stage of a two-stage disjunctive group testing 
algorithm under the condition that the vectors x are distributed with probabil- 
ities p\ x \l — p) v ~\ x \ where x G X(v) denotes the indices of the ones (defective 
items) in x. The bound above implies that the expected number of items which 
remain unresolved after application in parallel of u pools is not less than 

wEffW-pr-'a-* -v P . 
t=i ^ ' 
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3.3 Known Infinite Classes of Combinatorial Constructions 

Tonchev [25] straightforwardly gave a non-trivial construction method to obtain 
Steiner designs which have the additional property that the blocks have two sizes 
differing by one. 

Proposition 1 (Tonchev, 2008). Suppose that T> = (X,B) is a Steiner 
t-(v, k, 1) design that contains a Steiner (t — l)-(v, k, 1) subdesign T>' = (X, B'), 
where B' C B. Then, the blocks of T>' , each extended with one new point x (£ 
X , together with the blocks of T> that do not belong to T> , form a Steiner 
t-(y + 1, {k, k + 1}, 1) design. In particular, if there exists an RBIBD(v,k,l), 
then there exists a Steiner 2-(v + 1, {k, k + 1}, 1) design. 

Relying on resolvable BIBDs from afhne geometries and Kirkman triple sys- 
tems, Tonchev derived from the above result the following infinite classes: 

Theorem 3 (Tonchev, 2008). There exists 

• a Steiner 2-{q e + 1, {q, q + 1}, 1) design for any prime power q and any pos- 
itive integer e > 2, 

• a Steiner 2-(6a + 4, {3, 4}, 1) design for any positive integer a. 

Based on results on 2-resolvable Steiner quadruple systems by Baker [1] & 
Scmakov et al. [22] and by Teirlinck [23], Tonchev obtained this way also two 
infinite classes for t > 2. The third class had already been constructed earlier by 
Tonchev [24]. 

Theorem 4 (Tonchev, 1996 & 2008). There is 

• a Steiner 3-(2 2 f + I, {4, 5}, 1) design for any positive integer f > 2, 

• a Steiner 3-(2 ■ 7 e + 3, {4, 5}, 1) design for any positive integer e, 

• a Steiner 4-(4 e + 1, {5, 6}, 1) design for any positive integer e > 2. 

4 New Infinite Classes of Combinatorial Constructions 

We present several constructions of new infinite families of Steiner designs hav- 
ing the desired additional property that the blocks have two sizes differing by 
one. Our constructions involve, inter alia, resolvable BIBDs and cyclically resolv- 
able BIBDs. As a result, we obtain efficient two-stage disjunctive group testing 
algorithms suited for DNA library screening. 

4.1 CRBIBD-Constructions 

We obtain the following result: 

Theorem 5. Letp be a prime. Then there exists a Steiner 2-{pk + 1, {k, k + 1}, 1) 
design for the following cases: 
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(1) (k,p) = (3, 6a + 1) for any positive integer a, 

(2) (k,p) — (4, 12a + 1) for any odd positive integer a, 

(3a) (k,p) — (5,20a + 1) for any positive integer a such that p < 10 3 , and fur- 
thermore 

(3b) (k,p) — (5, 20a + 1) for any positive integer a satisfying the condition stated 
in (ii) in the proof 

(4) (k,p) = (7, 42a + 1) for any positive integer a satisfying the condition stated 
in (Hi) in the proof, 

(5) (k,p) — (9,p) for the values of p = 1 (mod 72) < 10 4 given in Table 2. 

Moreover, there exists a Steiner 2-(qk + 1, {k, k + 1}, 1) design for the fol- 
lowing cases: 

(6) (fc, g) for k = 3, 5, 7, or 9, and q is a product of primes of the form p = 1 
(mod k{k — 1)) as in the cases above, 

(7) (fc, g) = (4, q) and q is a product of primes of the form p = 12a + 1 with a 
odd. 

Proof. The constructions are based on the existence of a CRBIBD(pfc, k, 1) in 
conjunction with Proposition 1. We first assume that k is odd. Then the follow- 
ing infinite ((i)-(iii)) and finite ((iv)) families of radical difference families exist 
(cf. [3] and the references therein; [8]): 

(i) An RDF(p, 3, 1) exists for all primes p = 1 (mod 6). 

(ii) Let p = 20a + 1 be a prime, let 2 e be the largest power of 2 dividing a and 
let e be a 5-th primitive root of unity in 1 p . Then an RDF(p, 5,1) exists if 
and only if e + 1 is not a 2 e+1 -th power in Z p , or equivalently (11 + 5\/5)/2 
is not a 2 e+1 -th power in 7h v . 

(iii) Let p = 42a + 1 be a prime and let £ be a 7-th primitive root of unity in Z p . 
Then an RDF(p, 7, 1) exists if and only if there exists an integer / such that 
3^ divides a and e + 1, £ 2 + e + 1, e & re 3-^-th powers but not 3^ +1 -th 
powers in Z p . 

(iv) An RDF(p, 9, 1) exists for all primes p < 10 4 displayed in Table 2. 

Theorem 1 yields the respective CRBIBD(pfc, k, l)s. Moreover, in [12] a re- 
cursive construction is given that implies the existence of a CRBIBD(fcg, k, 1) 
whenever q is a product of primes of the form p = 1 (mod k(k — l j). In addition, 
a CRBIBD(5p, 5, 1) has been shown [4] to exist for any prime p = 1 (mod 20) 
< 10 3 . 

We now consider the case when k is even: In [19], a CRBIBD(4p, 4, 1) is 
constructed for any prime p = 20a + 1, where a is an odd positive integer. 
Furthermore, via the above recursive construction, a CRBIBD(4g, 4, 1) exists 
whenever q is a product of primes of the form p = 12a + 1 and a is odd. The 
result follows. □ 

Example 4- Values of p for which an RDF(p, k, 1) exists with k — 5, p < 10 3 , 
and k = 7 or 9, p < 10 4 are displayed in Table 2 (cf. [8]). For example, if we take 
an RDF(41, 5, 1), then we obtain a Steiner 2-(206, {5, 6}, 1) design. If we take an 
RDF(61, 5, 1), then we obtain a Steiner 2-(306, {5, 6}, 1) design. 



9 



Table 2. Existence of an RDF(p, fc, 1) with fc = 5, p < 10 3 , and fc = 7 or 9, 
p < 10 4 . 



k = 5 


41 61 241 281 


401 


421 601 641 661 


701 761 821 881 






k = 7 


337 421 463 883 


1723 


3067 3319 3823 3907 


4621 4957 5167 5419 


5881 


6133 8233 8527 8821 


9619 9787 9829 






k = 9 


73 1153 1873 2017 


6481 


7489 7561 



We remark that further parameters are given in [3] for RDF(p, k, l)s with k = 7 or 9 
and 10 4 <p< 10 5 . 



4.2 RBIBD-Constructions 

We establish the following result: 

Theorem 6. Let v be a positive integer. Then there exists a Steiner 
2-(v + 1, {fc, fc + 1}, 1) design for the following cases: 

(1) (fc, v) = (3, 6a + 3) for any positive integer a, 

(2) (fc, v) = (4, 12a + 4) for any positive integer a, 

(3) (k,v) = (5,20a + 5) for any positive integer a with the possible exceptions 
given in Table 3, 

(4) (fc, v) = (8, 56a + 8) for any positive integer a with the possible exceptions 
given in Table 3. 

Proof. The constructions are based on the existence of an RBIBD(w, fc, 1) in con- 
junction with Proposition 1. The following infinite series of resolvable balanced 
incomplete block designs are known (cf. [13,8] and the references therein): 

(i) When fc = 3 and 4, respectively, an RBIBD(u, fc, 1) exists for all positive 
integers v = fc (mod fc(fc — 1)). 

(ii) An RBIBD(u,5,l) exists for all positive integers v = 5 (mod 20) with the 
possible exceptions given in Table 3. 

(hi) An RBIBD(v, 8, 1) exists for all positive integers v = 8 (mod 56) with the 
possible exceptions given in Table 3. 

This proves the theorem. □ 

We remark that Case (1) has already been covered in Theorem 3. 
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Example 5. Choosing for example an RBIBD(65, 5, 1), we get a Steiner 
2-(66, {5, 6}, 1) design. If we choose an RBIBD(105, 5, 1), then we obtain a Steiner 
2-(106,{5,6},l) design. 



Table 3. Possible exceptions: An RBIBD(v, fc, 1) with fc = 5 or 8 is not known 
to exist for the following values of u = fc (mod fc(fc — 1)). 



k = 5 


45 


345 


465 


645 






k = 8 


176 


624 


736 


1128 


1240 


1296 1408 1464 


1520 


1576 


1744 


2136 


2416 


2640 2920 2976 


3256 


3312 


3424 


3760 


3872 


4264 4432 5216 


5720 


5776 


6224 


6280 


6448 


6896 6952 7008 


7456 


7512 


7792 


7848 


8016 


9752 10200 10704 


10760 10928 11040 11152 11376 11656 11712 11824 


11936 12216 12328 12496 12552 12720 12832 12888 


13000 13280 13616 13840 13896 14008 14176 14232 


21904 24480 











Theorem 7. If v and k are both powers of the same prime, then a Steiner 
2-(y + 1, {fc, fc + 1}, 1) exists if and only if (v — 1) = (mod (fc — 1)) and v = 
(mod fc). 

Proof. It has been shown in [14] that, for v and fc both powers of the same prime, 
the necessary conditions for the existence of an RBIBD(v, fc, A) are sufficient. 
Hence, the result follows via Proposition 1 when considering an RBIBD(t>, fc, 1). □ 

4.3 3-Design-Constructions 

Based on further results on 2-resolvable Steiner quadruple systems described in 
Teirlinck [23], we obtain this way also two infinite classes for t > 2. 

Theorem 8. There exists 

• a Steiner 3-(2 • 31 e + 3, {4, 5}, 1) design for any positive integer e, 

• a Steiner 3-(2 ■ 127 e + 3, {4, 5}, 1) design for any positive integer e. 
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5 Conclusion 

Efficient two-stage group testing algorithms that are particular suited for DNA 
library screening have been investigated in this paper. The main focus has been 
on novel combinatorial constructions in order to minimize the number of individ- 
ual tests at the second stage of a two-stage disjunctive testing algorithm. Several 
infinite classes of such combinatorial structures have been obtained. 
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